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Abstract 

As argued in the hterature, the heavy quark components in the collision hadrons may also 
provide sizable contributions to the hadronic production of the heavy flavored hadrons. Then, 
different to the fixed-flavor- number scheme (FFNS) as adopted in Ref. jld]. we study the hadronic 
production of and B* mesons under the general-mass variable-flavor-number scheme (GM- 

VFNS) so as to consistently deal with the double counting problem from the two mechanisms, i.e. 

- (*) - 

the extrinsic bottom mechanism via the subprocess g + b ^ Bs + s and the gluon-gluon fusion 
mechanism via the subprocess g + g^Bs +b + s. A comparative study of FFN and GM-VFN 
schemes are presented. At the TEVATRON, the difference between the two schemes is small, 
however such difference is obvious at the LHC, so the forthcoming experimental data on LHC 
shall provide a good chance to check which scheme is more appropriate to deal with the Bg meson 
production. 
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I. INTRODUCTION 



The properties of Bs meson, such as Bg — Bg oscillation, the mass difference Am 



and 



the charge-parity violation and etc., have been studied by CDF and DO collaborations 
Especially, the CERN LHC, with a high collision energy and a high luminosity, shall provide 
a good platform for studying the properties of Bg meson s], Q]. In the present paper, we 
shall focus our attention on its hadronic production. 

Principally, there are two approaches to deal with the heavy meson hadroproduction. One 
is the 'fragmentation approach', which is usually adaptable to higher transverse momentum 
region, e.g. Refs.[7-12| showed explicit examples for Be hadroproduction and 5-meson 
hadroproduction respectively. The other one is the so-called 'complete calculation approach'. 
Due to the fact that all the phase-space point can be conveniently retained within the 
'complete calculation approach', while the phase-space of the accompanied particles has 
been integrated out within the fragmentation approach so as to obtain the fragmentation 
function. Then the advantage of the complete approach is clear, i.e. it can retain the 
information of the accompanied quark jets, which sometimes is experimentally useful. Later 
on, we shall adopt the 'complete calculation approach' for our discussions. 

In the higher transverse momentum and smaller rapidity regions of a hadronic collider, 
which mainly corresponding to larger momentum fraction of the constitute quarks, it is 
found that the heavy quark components in proton or antiproton is always quite small in 
comparison to that of the light quarks or gluons. So in the literature, one usually does not 
take heavy quark components' contributions into consideration in most of the calculations 
for the heavy quarkonium production. For the purpose, the fixed-flavor-number scheme 



(FENS) is usually adopted 13Nl5|. Within FENS, the number of the active flavors in the 
initial hadron is fixed to be rij = 3, and then only light quarks/antiquarks and gluon should 
be considered in the initial state of the hard scattering subprocess. We have studied the 
hadronic production of the spin-singlet B,, and the spin-triplet B* within the 'complete 
calculation approach' by adopting FENS 161 ]. 

However, it has been argued that in certain cases the heavy quark components in the 
collision hadrons may also provide sizable contributions. For example, it has been shown that 
the mechanisms involving heavy quarks in the initial state can give sizable contributions to 
the hadronic production of (cc)-quarkonium [17|, (c6)-quarkonium [l8|, Scc-baryons jl9l. |20| . 



and etc.. Especially, it is found that in the lower pt region, the mechanisms from those heavy 
quarks in the initial state (hereafter refer to as "heavy quark mechanism" ) can even dominant 
over other mechanisms. The heavy-flavored quarks can be generated in two different ways 
inside the incident hadrons: firstly, it can be perturbatively generated by gluon splitting, 
and hence is usually named as 'extrinsic' component; and secondly, it can also be generated 
non-perturbatively and appears at or even below the e nerg y scale of heavy quark threshold 



that can be named as 'intrinsic' one according to Refs.|21l. |22j. For the present case of B'f^ 
production, the intrinsic 6-quark component in the hadron is quite small 23|, and it shall 
lead to negligible contribution for the production at both the hadronic colliders TEVATRON 
and LHC, so we shall not consider it at the present. While, as a comprehensive analysis of 
the hadronic production of B^*\ it would be interesting to study the extrinsic heavy quarks' 
contributions, which might be sizable, and also to make a comparison of the results with 
those obtained from FFNS. 



For such purpose, one has to deal with the production of 5^*^ through the heavy quark 
mechanism and the gluon-gluon fusion mechanism simultaneously. This is not an easy task, 
since we has to solve the double counting problem 2^: i.e. a full QCD evolved heavy quark 
distribution functions, according to the Altarelli-Parisi equations, includes all the terms 



proportional to In 



(yU stands for the factorization scale and mg is the heavy quark 
mass), while some of them also occur in the gluon-gluon fusion mechanism after doing the 
phase-space ir itegrat ion. We shall adopt the the general-mass variable-flavor-number scheme 
(GM-VFNS) 25I-I29I] to deal with the double counting problem. It should be noted that our 
present main purpose is to deal with the double counting problem between the heavy quark 
mechanism and the gluon-gluon fusion mechanism, and a full NLO calculation of the process 



is not available at the present due to its complexity, so our present imp 
VFNS is similar to the Aivazis-Collins-Olness-Tung (ACOT) approach 



ementation of GM- 



25j. 



The paper is organized as follows. In Sec. II, we present the calculation technology to 
estimate the i?^*-* hadropro duct ion under the GM-VFNS. In Sec. Ill, we present our numerical 
results and a comparison of the hadronic production of i?^*^ meson from the FFNS and the 
GM-VFNS shall also be presented. The final section is reserved for a summary. 
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II. CALCULATION TECHNOLOGY 



According to pQCD factorization theorem, the cross-section for the hadronic production 
of Bi*^ under the GM-VFNS can be formulated as 



+•••, (1) 

where the elhpsis stands for the terms in higher order and those terms with smaller 
contributions, such as the light quark and light anti-quark collision mechanism. The first 
term in Eq.([T]) shows the gluon-gluon fusion mechanism, which is dominant one under the 
FFNS. The second one represents the extrinsic 6-quark mechanism, in which the subtraction 
term is introduced to eliminate the double counting problem. The function F}j{x,fi) (with 
H = Hi or H2 and x = Xi or X2) stands for the distribution function of parton i in hadron H. 
da stands for the hadronic cross-section and da stands for the corresponding subprocesses. 
As conventional treatment, we have taken the renormalization scale fiR for the subprocess 
and the factorization scale for factorizing the PDFs and the hard subprocess to be the 
same, i.e. fiji = fip = fi. In the square brackets, the subtraction term for F^{x, /i) is defined 
as 

Flix, fi)suB = Ffj{x, /i) (g) /i) 

F]{n,^^)FfJC-.^^)-. (2) 

where the quark distribution Fg{x,fi) up to order is connected to the familiar g ^ bb 
splitting function Pg^ii ^ind its form can be written as 

Fhx, fx) = ^(1 -2x + 2x2) In 4- (3) 
^ 4tt mi 

b 

As for the gluon-gluon fusion mechanism, we need to deal with 36 Feynman diagrams at 



the leading order (a^), which has been dealt with in Ref. 16| by using the improved helicity 
amplitude approach as described in detail in Ref. 30|]. As for the extrinsic 6-quark mechanism 
for the hadronic B^*'> production, we need to consider the production via the subprocess 
9{P2) + b{pi) — )■ B'f\p^) + s{pi), whose typical Feynman diagrams for the extrinsic 6-quark 



P2 




Pi 




Pi X 

(6) 



FIG. 1. Typical Feynman diagrams for the heavy quark mechanism, i.e. hadroproduction of B, 
via the subprocesses: g{p2) + b{pi) — )• B^^\p^) + 5(^4). 



(*) 





FIG. 2. Graphical representation for the subtraction method within the GM-VFNS 25|. The 
symbol x on the internal quark line in the subtraction term indicates that it is close to the mass- 
shell and collinear to the gluon and hadron momentum. The combination of the first and the 
second terms are called as extrinsic mechanism. 

mechanism at the leading order (LO) are shown in Fig.([T]). In Eq.([T]), da r „(.) (xi, a;2, /x) = 

go ytSg 

c?(9"^j_^^(*)_(xi, X2, /i) stands for the usual 2-to-2 differential cross section, 



da 



{*) (xi,X2,/i) 



(27r)4|Mp 



n 



d^Pi 



^^{Pi-P2y- pipits (2^)'(2^> 



Pi -P2 



(4) 



where the initial-parton spin and color average and the final-state quantum number sum- 
mation are all attributed to |Mp. |Mp can be calculated by using the conventional way, 
only one needs to keep the heavy quark mass terms according to the GM-VFN scheme. 
For shortening the text, we do not put |M p here, whose explicit form can be found in the 
APPENDIX of Ref . [isj] by simply changing rric there to be the present case of m^. Next, 
to accomplish the phase space integration, we adopt two routines RAMBOS 3l|] and VE- 
GAS 32], which together with some reasonable transformations to make them run more 



33|, 



34 and GENXICC (35|,l36|. 



effectively can be found in the generators BCVEGPY [30, 

GM-VFNS is considered as the most complete framework for pQCD calculations being 
applicable in all regions of the phase space of hard processes. In particular, such scheme is 
by now standard in global analysis of PDFs to extract the most precise sets of partons to 
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be used for physics predictions at the LHC |37|. According to GM-VFNS, there is double 
counting of the gluon-gluon fusion mechanism and the extrinsic 6-quark mechanism. For 
the gluon-gluon fusion mechanism, in some of its Feynman diagrams, when the intermediate 
heavy quark line that is next to the incident gluon is nearly on shell and is coUinear to the 
incident gluon, then it will result in a factor of order distribution of a quark in a gluon, 
like Eq.(|3]). While such factor has already been included in the heavy quark distribution 
function, so these Feynman diagrams under such particular condition are double counted 
and must be subtracted. More explicitly, we draw Fig. ([2]) as an example to graphically 
illustrate this point, which is similar to the case of hadronic production of heavy quarks 



26(1 . The symbol (x) on the internal quark line in the subtraction term means that the 



heavy quark is on mass-shell and moving longitudinally. 

It should be noted that our present purpose is to deal with the double counting problem 
between the heavy quark mechanism and the gluon-gluon fusion mechanism, and a full NLO 
calculation of the process is not available at the present, so our present implementation 
of GM-VFNS is similar to the Aivazis-Collins-Olness-Tung (ACOT) approach [25], which 
sometimes is called as a LO generalized MS calculation 27|. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



In doing the numerical calculation, we take the decay constant and the radial wave 

0.209GeV and |i?(0)|2 = 0.497GeV^ 



function at the origin to be / 



B 



(*) 
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39| . The masses 



of h and s quarks are taken as mf, = 4.90GeV and nig = O.SOGeV. And to ensure the gauge 
invariance of the hard scattering amplitude, the mass of the bound state is taken to be the 



sum of the two constitute quark masses, i.e. M 



B 



(*) 



rrih + rris. The scale /i is set to be the 



transverse mass of the bound state, i.e. fi = Mt = ^jM^^^^^ + j^, where px is the transverse 
momentum of the bound state. Under the GM-VFN scheme, CTEQ6IIQ 40| is adopted for 



PDF, and to be consistent, the NLO as running above Aqqj^^^ = 0.326 GeV is adopted, i.e. 



An 



where /3o = 11 — 2nf/3 and /3i = 51 — 19ri//3. 



/3o 111(^2 /A^cd) 

As a comparison, for the FFN scheme, the LO as running is adopted. It is noted that for 
GM-VFNS, the value of changes with the energy scale. 




FIG. 3. pT-distributions of the extrinsic 6-quark mechanism for the hadronic production of Bg at 
LHC (Left) and TEVATRON (Right). The 'PURE-PDF' one means the 6-quark PDF is taken to 
be CTEQ6HQ, 'SUB.' one means the 6-quark PDFs is taken to be the subtraction term defined 
in Eq.dl]), 'COM.' one stands for the combination of the 'PURE-PDF' and 'SUB.' components as 
indicated by Eq.([T]). 



A. Hadroproduction of Bf' under the GM-VFNS 



Firstly, it would be interesting to show how double counting term is subtracted in GM- 
VFNS. We present the results of the extrinsic 6-quark mechanism for the hadronic production 
of Bs both at LHC and at TEVATRON in Fig. ([3]). It can be found that there exists a large 
cancelation between the contributions from the 'PURE-PDF' term (with 6-quark PDF taken 
to be CTEQ6HQ) and the subtraction term (with 6-quark PDF taken to be the subtraction 
term defined in Eq.(l2])), especially, in the large Pt regions. This shows that to obtain 
a reliable results under the GM-VFNS, we should take such double counting term into 
consideration, otherwise, the results maybe largely overestimated. 

Secondly, we present the total cross sections for the hadronic production of Bg and B* 
at LHC and TEVATRON, which are shown in Tab|T] and TabUl] respectively. Two typical 
cuts pTcut = 2.5 GeV and pxcut = 4 GeV for both LHC and TEVATRON, and \y\ < 1.5 
for LHC, \y\ < 0.6 for TEVATRON are adopted in the calculation. As a comparison, we 



also show the results under the FENS with PDF taken to be CTEQ6L [4l| in Tabjl] and 
Tab mi It should be noted that nj should be fixed to be 3 in the FENS and then to be 
consistent with the exact FENS, the PDFs for the initial partons should be taken the one 
like CTEQ5F3 [42,], which is generated by using the evolution kernels with effective flavor 



TABLE I. Cross section (in unit nb) for the hadronic production of at LHC {^/s = 14.0 TeV) 
and TEVATRON {^/s = 1.96 TeV) under GM-VFNS and FFNS, where the symbol g + b stands 
for the extrinsic 6-quark mechanism and etc. Two typical pT cuts are adopted. As for the rapidity 
cut, we take \y\ < 1.5 for LHC and \y\ < 0.6 for TEVATRON. 









GM-VFN 




FFN 




PTcut 


g + b 


9 + 9 


total 


9 + 9 


LHC 


2.5 GeV 


14.55 


44.43 


58.98 


83.06 




4.0 GeV 


4.41 


30.05 


34.46 


54.39 


TEVATRON 


2.5 GeV 


1.19 


3.28 


4.48 


4.50 




4.0 GeV 


0.32 


2.05 


2.37 


2.73 



TABLE II. Cross-section (in unit nb) for the hadronic production of B* at LHC {^/s = 14.0 TeV) 
and TEVATRON (^s = 1-96 TeV) under GM-VFNS and FFNS, where the symbol g + b stands 
for the extrinsic 5-quark mechanism and etc. Two typical pT cuts are adopted. As for the rapidity 
cut, we take \y\ < 1.5 for LHC and \y\ < 0.6 for TEVATRON. 









GM-VFN 




FFN 




PTcut 


g + b 


9 + 9 


total 


9 + 9 


LHC 


2.5 GeV 


32.94 


140.0 


173.9 


265.6 




4.0 GeV 


10.45 


92.24 


102.7 


167.9 


TEVATRON 


2.5 GeV 


2.54 


10.49 


13.03 


14.45 




4.0 GeV 


0.73 


6.33 


7.07 


8.47 



number ng// = 3. As argued in Refs. 43|, |4J], the uncertainties from different LO PDFs are 
small, and our numerically calculation shows that it only gives several percent difference by 
replacing CTEQ6L to CTEQ5F3 0, so as a conventional choice, we also adopt CTEQ6L as 
the typical PDF for FFNS. 

^ Under FFNS, by varying the flavor number with the energy scale, the value of as shall be decreased, but 
this is to a large degree compensated by a larger gluon distribution function (i.e. in small a;-region that 
is dominant for the production, F^(CTEQ6L) > Fj^(CTEQ5F3)), so as a whole, there is small difference 
by using CTEQ6L and CTEQ5F3. 
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FIG. 4. pt distributions for the hadronic production of Bg (Left) and B* (Right) at LHC under 
the GM-VFNS. The dash-dot and the sohd hues represent the gluon-gluon fusion mechanism and 
the extrinsic 6-quark mechanism respectively. Ah px distributions are drawn under \y\ < 1.5 and 
the PDF is taken as CTEQ6HQ. 




FIG. 5. Pt distributions for the hadronic production of Bs (Left) and B* (Right) at TEVATRON 
under the GM-VFNS. The dash-dot and the sohd hues represent the gluon-gluon fusion mechanism 
and the extrinsic 6-quark mechanism accordingly. All px distributions are drawn under \y\ < 0.6 
and the PDF is taken as CTEQ6HQ. 

It can be found from Tab|T] and TabHTl that the total cross sections of the extrinsic b- 
quark mechanism are comparable to those of the gluon-gluon fusion mechanism. The large 
cross section of the extrinsic 6-quark mechanism mainly comes from the small Pt region. In 
order to illustrate this point more clearly, we draw the transverse momentum distributions 
of B^/^ in Figs.dl]) and ([5]), which shows obviously that in a smaller pt region, the extrinsic 
6-quark mechanism is greater than that of the gluon-gluon fusion mechanism. However it 
drops down much more quickly with the increment of px, and then in a larger pt region. 




FIG. 6. pt distributions for the hadronic production of Bg (Left) and B* (Right) at LHC. The 
dashed and the dash-doted hues stand for the total (sum of the extrinsic 6-quark mechanism and 
the gluon-gluon fusion mechanism) results obtained under the GM-VFNS for rapidity cuts \y\ < 2.5 
and \y\ < 1.5 respectively. The solid and the dotted lines are for gluon-gluon fusion results obtained 
under the FFNS for rapidity cuts \y\ < 2.5 and \y\ < 1.5 respectively. 




FIG. 7. Pt distributions for the hadronic production of Bg (Left) and B* (Right) at TEVATRON. 
The dashed and the dash-doted lines stand for the total (sum of the extrinsic 6-quark mechanism 
and the gluon-gluon fusion mechanism) results obtained under the GM-VFNS for rapidity cuts 
\y\ < 1.5 and \y\ < 0.6 respectively. The solid and the dotted lines are for gluon-gluon fusion 
results obtained under the FFNS for rapidity cuts \y\ < 1.5 and \y\ < 0.6 respectively. 

the contribution from the extrinsic 6-quark mechanism is lower than that of the gkion-gluon 
fusion mechanism. 

Further more, it can be found that the gluon distribution of CTEQ6HQ is always smaller 
than that of CTEQ6L, especially in small x region, so the total cross section for the gluon- 
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gluon fusion under the GM-VFNS is smaller than that under the FFNS. Since x may reach to 
much smaller region at LHC than at TEVATRON, the difference between these two schemes 
is bigger at LHC than that at TEVATRON. Tab|T] and TabJITl show this point clearly. For 
example for the production of i?c, when pTcut = 4 GeV, it can be found that at LHC, the 
total cross section for the gluon-gluon fusion under the GM-VFNS is only 56% of that of 
FFNS; while at TEVATRON, such ratio raises up to 75%. This shows that when taking 
the extrinsic 6-quark mechanism into account for the GM-VFNS, one can shrink the gap 
between the GM-VFNS and the FFNS results to a certain degree. 

It can be found that in the large pt region, the extrinsic 6-quark mechanism is greatly 
suppressed due to the large cancelation from the subtraction term under the GM-VFN 
scheme, then the contribution from the gluon-gluon fusion mechanism shall be dominant 
in these higher px regions. In order to see the fact clearly, we present the pt distributions 
predicted by the GM-VFN and FFN schemes for the hadronic production of B^*^ at LHC 
and TEVATRON in Figs. ([6]) and (I7l) respectively. Figs. ([6]) and ([7]) shows that the main 
difference between the predictions by the GM-VFNS and the FFNS is only in small pt 
region {pt ^ 3.0 ~ 4.0 GeV). Furthermore, Fig. ([7]) shows that at TEVATRON, one can 
hardly distinguish the difference between FFNS and GM-VFNS, since a PTcut — 4:GeV is 
usually adopted at TEVATRON for hadronic productions. This shows that at TEVATRON, 
FFNS is enough to describe the data. While Fig. ([6]) shows that at LHC, such difference is 
amplified, so the forthcoming LHC experiment data may make a judge on whether we need 
to take the heavy quark component in proton into consideration (and hence the necessity of 
using GM-VFNS), since more small x and small pt events can be found/measured at LHC. 

B. A simple discussion on the uncertainties from 6-quark mass 

The uncertainties for the hadronic production of Bg and B* include the PDFs, the quark 
masses, the factorization scale and etc.. We have made a detailed discuss on such uncer- 
tainties ,„ Ref.Q ™der FFNS. here we shall o„.y eoncent.ate o„. atte„t,„„ o. the .qua* 

mass effect on the hadronic production, since the other uncertainties shall give similar be- 
havior under both GM-VFNS and FFNS @. For the purpose of discussing the uncertainties 

^ We have checked this conclusion numerically, and to short the paper, we do not present extra discussions 
on other uncertainties. 
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TABLE III. Cross-section (in unit nb) for the hadronic production of Bs[l^So] at LHC and TEVA- 
TRON under GM-VFNS and FFNS with nib varying within the region of (4.9 ± 0.1) GeV, where 
the symbol g + b stands for the extrinsic 6-quark mechanism and etc. As for the rapidity cut, we 
take \y\ < 1.5 for LHC and \y\ < 0.6 for TEVATRON. The upper, the center and the lower values 
are for m = 5.0, 4.9 and 4.8 GeV respectively. 







GM-VFN 


FFN 




PTcut 


g + b 


9 + 9 


total 


9 + 9 


LHC 


2.5 
4.0 


4.4llSJ^ 


44.43;|?2 
30.05;1:?2 


58.98;o;gg 

■^4 4fi-0.68 
"-"^•^"+0.77 


83.06;^!? 
54.39;!?^ 


Tevatron 


2.5 
4.0 


1 iq+0.12 
n OO+0.03 

"-'••3 -^-0.03 


q oo — 0.18 

2 05-°-°^ 


4-48+0.07 

r, Q7-0.06 

^•"3'+0.06 


4.5o;°:^^ 

r\ 'yq — 0.12 

^- '•3+0.13 



TABLE IV. Cross-section (in unit nb) for the hadronic production of i3*[l^S'i] at LHC and TEVA- 
TRON under GM-VFNS and FFNS with rrib varying within the region of (4.9 + 0.1) GeV, where 
the symbol g + b stands for the extrinsic 6-quark mechanism and etc. As for the rapidity cut, we 
take \y\ < 1.5 for LHC and \y\ < 0.6 for TEVATRON. The upper, the center and the lower values 
are for m = 5.0, 4.9 and 4.8 GeV respectively. 







GM-VFN 


FFN 






g + b 


9 + 9 


total 


9 + 9 


LHC 


2.5 
4.0 


32.94t|o™ 
in 4^+0-™ 

iU.40_Q g2 


140.0;6;l 
92.24;|9^ 


173.9;^ 
W2.7~li 


265.6;}2i 

167.9;!^ 


Tevatron 


2.5 
4.0 


^•"-'*-0.27 

n 70+0.06 
'•3-0.07 


io.49;°j^ 

fi S3-°-25 


lO.UO^Q 33 
'■'-"+0.20 


i4.45;°:i 
8.47;°i^ 



caused by rrib, we vary mi, within the region of (4.9 ± 0.1) GeV. The total cross section for 
the hadronic production of the scalar i?s[1^5'o] and the vector i?*[l'^5'i] at LHC and TEVA- 
TRON under GM-VFNS and FFNS are presented in TABUnland TABUVl For definiteness 
and focussing the uncertainties from the 6-quark mass alone, and we fix = 0.5 GeV, 
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fTT-Bs = f^s + '^b, f^'p = Ptb + ^^"-^ PDF is taken from CTEQ6HQ and CTEQ6L for 
GM-VFNS and PENS respectively. 

Erom TAB JIIII and TAB JIVt it is found that the total cross section for the extrinsic 
mechanism increases with the increment of mf,. Eor clarity, in the following, we shall take 
the case of Bg as an explicit examples to explain this point, and the case of B* is similar. 
When setting pTcut = 2.5 GeV, the total cross section for Bs shall be increased by 10% 
when ruf, increased by 0.1 GeV at both LHC and TEVATRON, which reduces to be 8% at 
LHC and 9% at TEVATRON for pxcut = 4.0 GeV. Inversely, since the allowed phase space 
becomes narrower with the increment of mi,, it is found that the total cross section for the 
gluon-gluon fusion mechanism decreases with the increment of rub. The total cross section of 
Bg for the gluon-gluon fusion mechanism shall be decreased by 5% (6%) when rrih increased 
by 0.1 GeV for prcut > 2.5 GeV at LHC ( at TEVATRON). As a combination, due to the 
different behavior of the extrinsic mechanism and the gluon-gluon fusion mechanism, the 
total cross sections under GM-VENS possess smaller uncertainties in comparison to that of 
EENS, which are 1.2% and ~ 1.5% at pxcut = 2.5 GeV when nib increased by 0.1 GeV at 
LHC and TEVATRON respectively. 



IV. SUMMARY 

As a compensation of our previous work [iG*] on the hadronic production of B^*^ meson, we 
have reanalyzed the production under the GM-VENS, and a comparison with the estimations 
of EENS and GM-VENS is also presented. Under the EENS, we only need to deal with the 



dominant gluon-gluon fusion mechanism as done in Ref. |16|. While under the GM-VENS, 
in addition to the gluon-gluon fusion mechanism, one also needs to consider the extrinsic 
heavy quark mechanism so as to make a more sound estimation. Eor the present case, it is 
found that the extrinsic 6-quark mechanism can be as important as the gluon-gluon fusion 
mechanism, especially in small pt region. This shows that if the hadronic experiments such 
as those at LHC and at TEVATRON can accumulate enough data and measure very low pt 
events then they can provide a good platform to check which scheme is more appropriate. 
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